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TECHNICAL NOTE No. 1549 

EFFECT OF YAW AT SUPERSONIC SPEEDS ON THEORETICAL 
AERODYNAMIC COEFFICIENTS OF THIN POINTED WINGS 
WITH SEVERAL TYPES OF TRAILING EDGE 
By W. E. Moeokel 


SUMMARY 

Expressions are obtained for the aerodynamic coefficients of 
thin pointed wings with three types of trailing edge. The deriva- 
tion is based on an equation previously derived for the pressure 
distribution on a thin yawed delta wing at supersonic speeds. For 
a semivertex angle of tan" 1 0.4, the effect of yaw on the lift 
coefficient, the wave-drag coefficient, the center of pressure, and 
the rolling-moment coefficient of a thin pointed wing is compared 
for three types of trailing edge. 

Although the lift coefficient of a pointed wing may be 
increased by sweeping the trailing edges back from the center line, 
the variation of center of pressure and rolling moment with yaw 
angle is thereby increased. The variation of center of pressure 
and rolling moment with yaw angle is decreased when the trailing 
edges are swept back from the outer tips of the wing. The chord - 
wise distance from the vertex of a symmetrical pointed wing to its 
center of pressure was found to be nearly independent of yaw angle 
for each of the types of trailing-edge sweepback investigated. 


INTRODUCTION 

Relations are derived in references 1 and 2 for the surface 
velocity potential of a thin-flat-plate yawed delta wing with 
leading edges Bwept behind the Mach oone. From these relations, 
expressions are obtained for the pressure distribution over the 
surfaces of the delta wing. For a fixed vertex angle, yaw angle, 
Mach number, and angle of attaok, the pressure coefficient was found 
to be constant along any radial line from the vertex to the trailing 
edge. When the pressure coefficient was integrated over the wing 
surface, explicit expressions were obtained for the lift coefficient 
and the center of pressure of the delta wing for’ any yaw angle for 
which the leading edges remain inside the Mach cone from the vertex. 
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The results of reference 1 are extended herein to include other 
types of pointed wing for whioh the pressure-coefficient equation is 
the same as for the delta wing. Equations are derived for the aero- 
dynamic coefficients of such wings when the trailing edge is swept 
back from the outer tips (wing i) , swept back from the center 
(wing II) , and swept back from both the center and the outer tips 
(wing III) . The lift coefficient, the wave-drag coefficient, the 
center of pressure, and the rolling-moment coefficient are evaluated 
for a semi vertex angle of tan -1 0.4 as functions of yaw angle and 
trailing-edge sweepback angle. 


SYMBOLS 


The following symbols are used in this report: 

Cj) wave -drag coefficient 

lift coefficient 

Cj rolling-moment coefficient 

Cp pressure coefficient 

E' complete elliptic integral of second kind, modulus V 1 " Gr 2 

F, G functions of 0;^, y, and 0 

M Mach number 

S surface area 

U free-stream velocity 

x, y coordinates parallel and normal to free-stream direction 

Xi, y^ coordinates parallel and normal to wing center line 

Xg, 72 coordinates normal and parallel to trailing edge 
x c maximum chordwise dimension 

x n normal distance from vertex to trailing edge 


x l> 71 
y 2 j 


> coordinates of center of pressure 
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a 

3 - 

¥ 

0 1 

e 2 

cp = tan 




angle of attack, radians 

1 / M 2 - 1 

angle of yaw, degrees 

semiangle of vertex of pointed wing, degrees 
angle between trailing edge and y^-axis, degrees 
tan ©2 \ 


■^a 
A b 
®a 
®b 
a = 


tan tan ©gy 

= tan ©i + tan ©g - (l - tan ©2 tan ©g) tan(©g - cp ) 

= tan ©^ + tan 0 g - (l - tan ©2 tan ©g) tan(©g +cp) 

= tan ©^ - tan ©g + (1 + tan ©2 tan 0 g) tan(©g - cp ) 

= tan ©2 - tan © 2 + (l + tan ©2 tan ©g) tan(©g +cp) 

tan ©n 


a' = 


tan ©2 + tan f 
tan ©2 


b = 


b' = 


tan ©2 - tan y 

1 

tan ©2 + tan y 

-1 


tan ©2 - tan ¥ 


c = a - b tan © 2 
c ' = a ’ - b ' tan ©g 
e = a tan ©g + b 
e 1 =s a 1 tan ©g + b 1 

_ 2aFK _ 

~ r 2 - r2 

12, Ig, I3 definite integrals used in evaluating C j, and center of 
pressure 
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K = tan(0-L +?) tan (9^ - ?) 

Kl - — i- 

^ -ee ’ 

Ko = a/ iJL*” 8 1 

1/ 1 + tan 0 X 


tan ©2 


tan 02 


K 3 = f\j(o + e tan 02 ) (o' + e' tan 02 ) 

^ ce 1 + c ’ e 
*4 = g 


K 5 = 3 K 4 K 1 " 

Kg — K 4 Kg + cc 1 

r * integration parameter, Vz/'X-Z 

ri lower integration limit 

r 2 upper integration limit 

Subscripts: 

a, b, c, d segments of wing 

¥ evaluated for jiositive values of 2' 

-2f evaluated for negative values of ? 


ANALYSIS ‘ 

In the notation of this paper, the equation for the pressure 
coefficient Cp of a thin delta wing at supersonic speeds, as 
derived in reference 1, Is 



where K, a, b, a', and b' are functions of the semivertex angle 0-j. 
and the yaw angle ¥> and yj. and x^ are coordinates fixed with 


&fr8 
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respect to the wing. (See fig. 1.) The method used in referenoe 1 
to derive equation (l) was approximate. The approximation, however, 
affected only the values of the factor F. An exact expression 
for F has been obtained in reference 2. In the notation used herein, 
this expression is 


F 



2 Q. 

) + tan (©! - ¥ )] 


3 [tan (©]_ + 


( 2 ) 


where 


G = 


l+p 2 tan(e 1 +¥)tan(e 1 -'S)- f\j [l-3 2 tan 2 (©-,_+¥)] [l-p 8 tan 2 (e 1 -?5 


3 [tan (&i + 5 ) + tan (Sq - ? J] 


Values of F for several values of ©q and for two values of 3 
are shown in figure 2. A plot of the variation of pressure coeffi- 
cient with yaw angle tjr over a wing with semi vertex angle of tan ” 1 0.4 
is shown in figure 3, which is taken from figure 6 of reference 1. 

(The use of equation (2) in place of the approximate values of F had 
no noticeable effeot on the values shown in fig. 3.) 

Because Cp is constant along a given line from the vertex, 
equation (l) is also valid for the types of wing shown in figure 4 
provided that the leading edges fall behind the Mach waves from the 
vertex and the trailing edges fall ahead of the Mach waves from the 
foremost points of the trailing edges. These conditions may be 
written as 


©1 + ijr —cot ” 1 3 (3) 

© 2 + ? ^-90° - cot " 1 3 (3a) 

Equations (3) and (3a) imply that 

e 1 + ©2 +25 ^90° (3b) 

In order to simplify the integrations required to determine 
lift coefficient and center of pressure, it is convenient to obtain 
equation (l) in terms of coordinate axes parallel and perpendicular 
to the trailing edge (fig. 4) . With this coordinate system, the 
integration need be performed for only one variable and the xg 
coordinate of the center of pressure is oonstant (xg = 2x a /3) . 
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Xq_ « Xg cos ©2 + ^2 s ^ n ®2 


yi 

y± 

x i 


72 cos ©2 " *2 8in ®2 

VZ 

tan ©o 

x 2 <= 

yi 

i + — tan © 2 


(4) 

(4a) 

(4b) 


When equation (4b) Is substituted into equation (l), the following 
expression is obtained for the pressure coefficient: 



where r = y2/ x 2* Tbe lift coefficient and the oenter of pressure 
are then obtainable from the following equations: 
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Because the definite integrals in equations (6) and (7) occur for 
each of the wings, these equations are written in the following form.: 

C L = H [(o +c*)li +(©+©') I 2 J ( 6 a) 

— = l"(c + c*) Ig + (e + e*) I 3 I (7a) 

x 2 C L L J 



in order to obtain the lift coefficient and the center of pres- 
sure for the wings of figure 4, it remains only to evaluate I p, Ig, 
and I3 for the proper integration limits for each wing. The 
expressions for the quantities to he used in equations (6a) and (7a) 
sure summarized for each wing in table I. The limits r2 and rp 
may be found in the expressions given for H/a. 

For wings I and II, the expressions in table I are the same for 
either half of the wing. The lift coefficient for the unsheided half 
in figure 4(a) or 4(b) is obtained when ? is replaoed by -? in 
calculating the values of c + o' and e + e' in equation (6a). 

The lift coefficient for the entire wings is one-half the sum of the 
lift coefficients of the two halves; that Is, 
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°L = \ [( c L,a)«[» + ( 12 ) 

where (C^ a)_f is the lift coefficient obtained from equation (6a) 
when - $ is substituted for $ in the terms c + c ' and e + e ' . 
Thus by carrying out the integrations over half of the wing surface 
and calculating the lift coefficient for negative as well as positive 
values of ? , the lift coefficient of the entire wing may be obtained. 


Coordinates of the center of pressure for the entire wing are 
similarly obtained from values computed for only one-half of the 
wing by reconverting the coordinates obtained from equation (7a) to 
the x-^, yq coordinate system. Then, for the entire wing 
(fig§. 4(a) and 4(b)), 


- _ (°L,a yi, a )_? 

^1 

(°L,a) ? + ( c L,a)-$ 


(13) 


5 = (C L,a x l,a>g + ^ C L,a x l, a *-g 

( c L,a)f + ^ c L,a)-f 


(13a) 


For the wing of figure 4(c), the expressions in table I 
pertain to the segments a and b. The lift coefficients for each 
half of the wing are the mean of the coefficients of the two 
segments weighted according to the areas of the segments; that is 


^a+c 


c L,b+d = 


s a(°L,a)5 + Sb^I^b)-? 
®a + ®b 

s a( c L,a)-$ + s b( c L,b) , P 
s a + Sb 


(14) 


(14a) 


The coordinates of the centers of pressure for each half of the 
wing are: 


s a( c L,a yi,a)? " ®b^ C L,b yi,b)-? 


yi,a+c 


s a(°L,a)? + ®b^L,b^- V 


(15) 
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_ s a(°L,a x l,a)? + M^b x 1,1)^-¥ 

x l,a+c 

s a(°L,a)? + ®b( c L,b)-¥ 


- s a(°L,a yi a >-¥ + ®b^L,b y l,b)f 

yi,T)4d - 

Sa(°L,a)-? + Sb^L^b)? 


(15a) 


(16) 


x l,b+d 


ga(PL,a x l,a)-g + gb(^b x l,b) f 
s a( c L,a)-¥+ ^(^b)? 


(16a) 


where 


S » - [ t<m (fl 2 



- cp ) - tan (©2 " ®i)J 

(17) 

+ cp) - tan © 2 J 

(17a) 


In order to obtain the lift coefficient and the center of pres- 
sure for the entire wing, the values obtained for each half are 
combined as in equations (12) and (13). For 0g ** equations 

for the lift coefficient and for the center of pressure for all three 
wings reduce to those obtained for the delta wing in reference 1} 
that is, 

C = 2rtog tan 9 ± (l8 ) 

^1 - tan 2 tan 2 ¥ 


71 tan W 



(19) 

(19a) 


RESULTS 

The results of computations for pointed wings having a semi- 
vertex angle of ©q = tan -1 0.4 » 21.8° are presented in fig- 
ures 5 to 7. 
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The lift-slope curves (which are also the drag-slope curves 
divided by a, if leading-edge suction is neglected, because 
C D /C L = a for the flat-plate wing) are plotted in figure 5 as 
functions of yaw angle and trailing-edge sweepback angle for a 
Mach number of / \f2. Results are presented for half of each wing and 
for the entire wing. The variation of lift coefficient with yaw 
angle is similar for each of the three wings considered. As the 
trailing-edge sweepback angle 62 is increased, the lift coeffi- 
cient decreases for wings I and III, but increases for wing II. 
These results are to be expeoted because, as shown in figure 3, 
the pressure distribution is assymptotic to the outer edges of 
the wing. Hence, removing portions of the wing from the outer 
edges should lower the lift coefficient. The variation of lift 
coefficient with 9q is slightly greater for wing III than for 
wing I.‘ 

The variation of the center of pressure with yaw angle and 
trailing-edge sweepback angle is shown in figure 6 for the wings 
of figure 4. Because the factor F, which is the only quantity 
that varies with Mach number, cancels in the calculations of 
center of pressure, these results are valid for all Mach numbers 
so long as the conditions of equations (3) and (3a) Eire fulfilled. 
The variation of center of pressure with yaw angle is shown for 
one-half of each wing as well as for the entire wing. The results 
for one-half of each wing, of course, have significEtnce only as an 
aid in computing the results for the entire wing. For halves of 
wings I and II, the center of pressure moves along lines parallel 
to the trailing edge as the yaw angle is varied. For half of 
wing III, the xp coordinate of the center of pressure Is only 
slightly dependent on the yaw angle. 

For the complete wings I, II, and IU, the X]_ coordinate of 
the center of pressure remains constant as ? Is varied. For 
fixed values of ?, the center of pressure for each of the wings 
moves along straight lines as Sg * s varied. In particular, for 
wing I (fig. 6(b)) the lines for fixed values of ? are inclined 
at the angle to the center line and appear to intersect at a 
common point whose coordinates are y^ = 0, x^ = x c /3 . Hence, the 
approximate equation for the variation of center of pressure with 
yaw angle for this wing is 


n 

*c 


tan * (S . i) 


( 20 ) 
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For wings II and III, the center of pressure for fixed values of f 
moves along lines that appear to have a common intersection point 
at the coordinates y]_ = 0, xj_ = 1.43 x 0 and y^ *= 0 , % = 0.473 !<., 
respectively. For these two wings, however, the lines are not 
inclined at the angle ? with respect to the xq/x c axis. The 
reasons for common intersection points are not obvious from the 
equations for the center of pressure and these intersections may "be 
only apparent, that is, only true within the accuracy of the plots. 


The intersections of the constant lines for halves of 

wings I and II at the coordinates yq = 0, xq = 2xq/Z and 
yq = 0.4xq, Xq = 2x c /3, respectively, also have no readily 
discernible reason although a common Intersection point is plausible 
because each of these lines mist be parallel to the trailing edge. 
(The constant ©3 lines are parallel to the trailing edges because 
the center of pressure for eaoh elementary area of the half -wings 
lies 2/3 of the distance from the vertex to the trailing edge.) 

A comparison of figures 6 (b), 6 (d), and 6 (f) shows that the 
center of pressure varies considerably more with yaw angle for 
wing II than for wings I and III. This result is to be expected 
because the pressure is highest along the edges of the wings and 
least near the center. (See fig. 3.) Removing portions along the 
outer edges should therefore reduce the center-of -pressure travel. 
The lift coefficient, of course, was also less for wings I and TH. 


The variation of the rolling-moment coefficient with yaw angle 
and trailing-edge sweepback angle is shown for wings I and II in 
figure 7. This coefficient is defined by the equation 




c L 


yi 


and is a measure of the rolling moment about the center line of the 
entire wing. For wing I the variation of Cj with yaw angle is 
less than for wing II. The variation of rolling-moment coefficient 
with yaw angle and ©2 vas almost the same for wing III as for 
wing I because the lift coefficient is slightly smaller and the 
center-of -pressure travel slightly greater for wing III than for 
wing I. 

Because the xq coordinate of the center of pressure remains 
constant for each of the wings for a given © 3 > no change in 
pitohing moment with yaw angle need be expected under steady-state 
conditions for any of the wings considered if the center of gravity 
is on the center-of -pressure line obtained for 5 = 0 . 
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CONCLUSIONS 

From calculations of the variation of aerodynamic coefficients 
of pointed wings with yaw angle and trailing-edge sweepback angle, 
the following conclusions may be drawn: 

1. The lift coefficient of a pointed wing may be increased by 
sweeping the trailing edges back from the center line. 

2. The chordwise distance from the vertex of a symmetrical 
pointed wing to the center of pressure is nearly independent of 
yaw angle for any of the types of trailing-edge sweepback investi- 
gated. 


3. The variation of center of pressure with yaw angle is 
decreased when the trailing edges are swept back from the outer 
tips and is increased when the trailing edges are swept back from 
the center line. 

4. The variation of rolling-moment coefficient with yaw angle 
for a pointed wing may be reduced by sweeping the trailing edges 
back from the outer tips. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 30, 1947. 
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TABLE I« - VALUES REQUIRED TO OBTAIN LIFT COEFFICIENT AND CENTER OF PRESSURE FOR WINGS OP FIGURE 4 


Wing I, segment a | Wing II, segment a 


2FK 


2K^ tan*l Kg 


tan T§2 +&x) - tan 0g 


2Ki(f - tan' 1 Kg) 


Wing III, segment a 


2FK 


tan (eg-*) -tan (Bg-Si 


Wing III, segment b 


2FK 


tan (e 2 +*) - tan 9g 


“i ^WS 




I 2 -K* (tt 3 - K4 1^ Kl 2 (K 3 *K 4 I ^ 


tan £ '9x-tan fa T 




tarr^-tarr 


_*q a J~| A b B b 


■ - S 3- K 4 1 


cos &2 


cos (8g-»9i) 


cos 0i 


cos (0g-*) 


cos (e 2 +*) 
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(b) Wing II; trailing edges swept back from center. 

Figure 4. Continued. Geometric parameters for pointed wings with swept' 
back trailing edges. {Equations in table I apply to shaded areas.) 







Gj/a and Cjj/«^ 
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Figure 5- - Lift and drag coefficients as function of yaw angle and 
t ral I I ng-edge sweepback angle 9 \ “ tan - * 0.4; M - */57 
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gure 5- - Continued- Lift and drag coefficients as function of 
angle and trai 1 ing-edge sweepback angle 9 9 ^ = tan - 0.4; M 


yaw 

- -f2. 
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'g —24 -16 -8 0 8 16 24 

« Yaw angle, f , deg 



Yaw angle , ± ¥ , deg 
(f) Entire wing III. 


Figure 5. - Concluded. Lift and drag coefficients as function of y aw 
angle and t ra i I i ng-edge sweepback angle 8 2 * 8 | « tan - ^ 0.4; M — V27 
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(b) Entire wing I. 

Figure 6. - Variation of center of pressure With yaw angJe and trai ling- 
edge sweepback angle- 9 { = tan“ * 0-4- 
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Figure 6- - Continued. Variation of center of pressure with yaw angle 
and t ra i I Ing— edge sweepback a ngj e. 9 j — tan - 0.4. 
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Figure 6. 
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(f) Entire wing III. 

- Concluded. Variation of center of pressure with yaw angle 
and trai 1 ing-edge sweepback angle. 8 j * tan - ^ 0.4. 











